Zero-temperature dynamics of solid ^He 
from quantum Monte Carlo simulations 
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The lattice dynamics of hep crystalline *He is studied at zero temperature and for two different 
densities (near and far from melting), using a ground-state path- integral quantum Monte Carlo 
technique. The complete phonon dispersion is obtained, with particular attention to the separation 
of optic and acoustic branches and to the identification of transverse modes. Our study also sheds 
light onto the residual coherence affecting quasiparticle excitations in the intermediate momentum 
region, in-between the phonon and nearly-free particle regimes. 



The lattice dynamics of solid ^He has long been con- 
sidered a major challenge to ab initio calculations, due 
to the strong anharmonicity of this highly quantum solid. 
Many calculations have been performed within the Self 
Consistent Phonon (SCP) approximationfi*^ which how- 
ever may be rather unsatisfactory, due to the magnitude 
of anharmonic effects. A significant improvement has 
been made possible by the application of a variational 
quantum Monte Carlo (QMC) approach, based on the 
shadow wave function formalism, to bcc ^He^ and to hep 
^Hgj^ The variational nature of this approach, however, 
makes it not fully suitable at high energy where optical 
or zone-boundary longitudinal excitations exhibit broad 
multiphonon features. In this spectral regime, the cal- 
culation of the full dynamic structure factor is therefore 
in order. QMC techniques based on path integrals allow 
quite naturally for the calculation of imaginary-time cor- 
relation functions from which various spectral functions, 
such as the dynamic structure factor, can be obtained 
upon analytical continuation. This technique has been 
successfully demonstrated for supcrfluid "^Hc^^ as well 
as for the bee crystalline phases of ''He — and ^Hci^ In 
the latter studies the spectrum of the transverse excita- 
tions has also been obtained, albeit in the one-phonon 
approximation onlyi^ 

In this work we present an extensive study of the dy- 
namical properties of hep ^He at zero temperature, per- 
formed by estimating the dynamic structure factor from 
ground-state path-integral simulationS] ^i^°'^^ This tech- 
nique allows us to parallel to some extent the procedure 
followed experimentally to map phonon dispersions from 
the measured neutron scattering. In the long wave-length 
region — well approximated by a phonon picture of the 
collective density excitations — we thus obtain longitudi- 
nal as well as transverse modes for both acoustic and 
optical branches. For higher wave- vectors we analyse the 
dynamic structure factor in terms of corrections to the 
so-called impulse approximation^^ finding a coherent re- 
sponse which is peculiar of both superfluid and solid he- 
lium. 

In Sec. Uwe give an introductory account of the phonon 
theory of long wave-length excitations in solids. In Sec. 
|TT]thc reader is provided with an outline of the numerical 



methods adopted in this work. In Sec. IIIII we report 
on the analysis of our QMC results both in the phonon 
regime and in the intermediate momentum region. Sec. 
IIVI is finally devoted to a few concluding remarks. 



I. LATTICE DYNAMICS 

A. Long wave-lengths 

The long wave-length lattice dynamics of a solid is fully 
characterized by its dynamic structure factor, which is 
the space-time Fourier transform of the density-density 
correlation function. In real time and reciprocal space, 
the autocorrelation function of the density operator 
reads: 
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where the brackets indicate equilibrium (ground-state or 
thermal) expectation values. In a weakly anharmonic 
system, it is convenient to expand S'(Q,t) into a sum of 
terms involving one-phonon processes, two-phonon scat- 
tering, interference processes and so on;^ 



SiQ, t) = 5i(Q, t) + 52(Q, t) + SiM t) 
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The physical meaning of such an expansion is best ap- 
preciated by introducing the atomic displacements from 
the equilibrium lattice sites, {R;}: u;(t) = ri{t) — R;. In 
terms of the u's and the R's, the one-phonon contribu- 
tion to the dynamic structure factor of a simple Bravais 
lattice reads!^ 
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where e~^^ is the Debye- Waller factor. For a harmonic 
crystal — to which only, strictly speaking, the phonon 
language applies — we consider the vibrational frequency 
tjj(q) and polarization vector e(q|j) of the j-th phonon 
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branch at wave vector q in the first Brillouin Zone (FBZ). 
In terms of these quantities, the one-phonon contribution 
reads )^ 
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where Q = q+G, G being a reciprocal-lattice vector, and 
.9^(Qb) IQ • e(QI.?)l^ is the so-called inelastic structure 
factor that filters out transverse vibrations. In the case 
of a non-Bravais lattice, such as the hep phase of Helium, 
the form of the inelastic structure factor is slightly more 
complicated^^ 



At intermediate wave-lengths both the phonon and a 
purely impulsive, particle-hole, picture of density exci- 
tations break down. In spite of the attention paid by 
both experimentalists and theorists to this peculiar in- 
termediate regime, both in the superflui d^^i^^ and in 
the soli d^^i^° phases, it turns out that the neglect of 
interaction-induced coherence effects make previous the- 
oretical studies not totally satisfactoryii^ 

At small wave-length, the solid behaves like a collec- 
tion of almost non-interacting atoms and the intermedi- 
ate scattering function can be approximated by its inco- 
herent partfi^ i.e. 
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where the d's are the positions of the atomic basis. 

In a perfectly harmonic solid the Fourier transform 
of Eq. (U)), 5i//(Q,ti;) is merely a sum of Dirac delta 
functions centered at the phonon frequencies. In a real 
solid, things are more complicated: anharmonic interac- 
tions broaden the one-phonon peaks and give rise to non- 
vanishing multi-phonon and interference contributions to 
the dynamic structure factor (Eq. [5]). When anharmonic 
effects are not too large, one-phonon excitations can still 
be long-lived — thus providing a reasonable description 
of the dynamics — and it is thus well justified to iden- 
tify the positions of the finite- width peaks of S'i(Q,w) 
with phonon frequencies. From an experimental point 
of view, phonon frequencies are generally extracted from 
the peaks of the full dynamic structure factor S'(Q,tt'). 
The cross section of inelastic neutron or X-ray scatter- 
ing is in fact proportional to 5(Q,a;)^ and no direct 
access is possible to its one-phonon component. The lat- 
ter dominates the cross section only at small transferred 
momentum, whereas multi-phonon contributions cannot 
in general be neglected when pursuing a comparison be- 
tween calculated and measured phonon dispersions. 



B. Shorter wave-lengths 

The very concept of phonon^ which lies at the basis 
of the theory of lattice dynamics sketched above, is most 
appropriate to describe the low-lying portion of the spec- 
trum of solid ^He, probed by inelastic neutron or X-ray 
scattering at long wave-lengths. In the opposite limit 
of short wave-lengths, the scattering process can be pic- 
tured as the creation of particle- hole pairs, resulting from 
the high momentum transferred to the crystal from the 
incoming particle beami^ The kinematics of the struck 
particles in the higher-energy states is clearly affected 
by the distribution of allowed atomic momenta, 7t,(Q), 
and neutron spectroscopy at large momentum transfer 
has in fact proven useful to probe off-diagonal long-range 
order, both in superfiuidiS. and, more recently, in solid 
Hefium-i^ 
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which amounts to neglecting the interference terms in- 
volving different atoms. For a crystal, the incoherent 
part can be expressed in terms of the recoil frequency 
ujii = 2m^^ ^^"^ '^^ phonon density of states ^(w), 
leading to 
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Such an expression has been used by Glydei^ to com- 
pute the incoherent response of hoc ^He. By its very 
nature, the incoherent approximation is only reliable at 

very high wave- vector roughly larger than 20 A .In 
order to account for the leading coherence effects on the 
short wave-length dynamics of an extended system, it is 
convenient to consider a cumulant expansion of the in- 
termediate scattering functioufi^ 
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where S'(Q) is the static structure factor and are the 
cumulants of the distribution S'(Q,tt' — uir). Retaining 
the leading contribution to such an expansion yields the 
so-called Impulse Approximation^ according to which the 
Fourier transform of the dynamic structure factor con- 
sists of a main Gaussian component centered at the recoil 
frequency loh, plus additive corrections)^ 
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where the first terms of the expansion read 
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with oj^ = ujii/ S{Q,) and tj^ = (cj — ti-"^)^ //^2 



II. NUMERICAL METHODS 

The dynamical properties of bosonic systems are con- 
veniently simulated in imaginary time, r = it^ using path 
integral QMC methods, at both finite^ and zero^ temper- 
ature. Specializing to the zero temperature case, a dis- 
cretized path integral expression for the imaginary-time 
propagator can be used to project out the exact ground 
state ^'o from a positive trial wave function <f>o, according 
to 



The trial function is of the standard McMillan- 
Nosanow form, 



MR) 



exp 



L k<l 



(12) 



where {ri,...,rjv} = R are the coordinates of the N 
atoms, and s; is the i-th site of the hep lattice. The three 
parameters a, b, and c are optimized by minimizing the 
variational energy and their numerical values for the den- 
sities considered in this paper arc shown in Table |T1 The 
Gaussian localization terms in the trial function explicitly 
break Bose symmetry at the variational level. Nonethe- 
less, we believe that the lack of permutation simmetry 
cannot affect the determination of the density excitation 
energies due to the small frequency of particle exchanges 
in the crystal. This point can be further elucidated notic- 
ing that the presence of a vacancy in the simulation box 
greatly enhances the number of exchanges Therefore, 
if indistinguishability were important, one would expect 
a change in the frequency upon doping with vacancies. 
However, no such effect was found in in a variational 
calculation^ where Bose symmetry was taken into ac- 
count. 
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thus mapping the imaginary-time evolution, from which 
ground-state expectation values can be obtained, onto 
a classical system whose fundamental variables are open 
quantum paths (or reptiles in the parlance of Refs. and 
[22! ) . Full details of the formalism can be found in Ref. [13 
and need not be repeated here; we only stress that unbi- 
ased ground-state expectation values are obtained when 
the projection time (3 is large enough and the step e of the 
time discretization is small enough. We simulate samples 
of ^He atoms interacting through the Aziz^ pair poten- 
tial, placed in a cuboid cell accommodating an hep lat- 
tice. The number N of particles is either 180 or 360, the 
latter corresponding to a cell with double extension in 
the TA direction. Altough all the reported results refer 
to the larger system, we have found a full agreement be- 
tween the relevant observables computed on the common 
set of wave vectors shared by the smaller and the larger 
simulation cells. For systems of this size, we find it more 
efficient to use the bisection algorithmiSiii rather than 
the reptation algorithm^ for sampling the path space. 
We adopt the so-called primitive approximation for the 
imaginary-time propagator^ which requires a small time 
step e = 10^'^ inverse K for accurate results, and we set 
the projection time to /3 = 0.3. The correlation func- 
tions of the density operators are calculated for values of 
r up to Tmax = 0.25, implying that our paths span an 
imaginary time of 2(3 + Tmax = 0.85. 



0.028 156.0 5.74 0.565 
0.037 139.4 5.75 0.993 



Table I: Optimal values of the variational parameters appear- 
ing in the McMillan-Nosanow wave-function. 



In order to obtain ^(Q,^) from imaginary-time corre- 
lations, an inverse Laplace transform must be performed, 
for which we use the maximum entropy methodi^ Al- 
though the reconstructed spectra are typically much too 
broad, this procedure gives good results, at least for the 
position of the peaks, when a single sharp feature ex- 
hausts most of the spectral weight. As a typical exam- 
ple, we show in Figure [T] two spectra at different wave- 
vectors for the longitudinal acoustic branch, calculated 

° —3 

at the melting density pi = 0.028 [A ]. Despite the 
fact that much of the width of the peaks is an artifact 
of the numerical inversion of the Laplace transform, it is 
nonetheless plausible that the broadening of the spectra 
shown in Figure [T] reflects stronger multiphonon effects 
at higher wave-vectors. In the following we refer to the 
position of the peaks of the reconstructed spectrum as to 
phonon energies. The reported error bar is the statistical 
uncertainty of the peak position, as estimated with the 
jackknifc resampling method 
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Figure 1: (Color online) Example of two reconstructed spectra 
at different wave- vectors along the FA direction at the melting 
density pi. 



III. RESULTS 

In this Section we present the results of our QMC sim- 
ulations for both the phonon dispersion energies and the 
higher wave- vectors response of the solid. 



A. Long wave-length excitations 

Longitudinal modes 

The excitation energies of longitudinal vibrations can 
be straightforwardly obtained from the dynamic struc- 
ture factor 5(Q,a;), when available. In a lattice with 
a basis, such as hep ^Hc, multiple branches (acoustic 
and optic) exist at each point of the FBZ. Although 
for a generic wave-vector all the branches contribute to 
5(Q,tt'), it often happens that — because of the explicit 
dependence of the inelastic structure factor [5] on both the 
wave- vector and the branch index — along high-symmetry 
directions different branches dominate (and in practice 
are only visible) at different values of the wave- vector. As 
a consequence, there are regions of the reciprocal space 
in which the acoustic modes dominate while the optical 
modes are suppressed and vice-versa. To figure out the 
relative weights of the branches, the inelastic structure 
factor can be calculated in a number of (approximate) 
ways, as in Ref. [s^ for beryllium, another hep solid. By 
virtue of the strongly geometrical nature of 5^(Q|j), it is 
sufficient to look at one of these approximate calculations 
performed for the hep geometry to realize that, with few 
exceptions, the relative weights do generally suppress one 
mode and privilege the other. Our results substantially 
confirm this picture, the calculated spectral functions be- 
ing generally dominated by a single peak. Reconstructing 
a complete picture of the phonon dispersions thus require 
sampling the dynamic structure factor outside the FBZ. 



We have calculated the phonon energies at the melting 
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density pi = 0.028 [A ] — in a regime of strong quan- 
tum fluctuations signaled by a considerable Lindcmann's 
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ratio — and at the density p2 = 0.037 [A ], where the 
quantum fluctuations are less pronounced. Results are 
shown in Figures [2] and [3] using an extended-zone scheme 
reminiscent of the way the optic and acoustic modes are 
measured in the lab. The phonon energies extracted by 
5(Q, uj) are compared to experimental data. The phonon 
energies resulting from an analysis of the one-phonon con- 
tribution to the dynamic structure factor are also shown 
for comparison. 

The main findings that emerge from the calculations 
of the longitudinal modes are the following: 

1. The overall agreement between the calculated and 
measured peak energies is good. The estimated er- 
rors come from the intrinsic width of the peaks of 
the dynamic structure factor, which is larger for op- 
tic than for acoustic phonon. This feature is present 
in both the theoretical and experimental spectra, 
although in the former the width is enhanced by the 
numerical difficulties in performing inverse Laplace 
transforms. 

2. The discrepancy between the frequencies estimated 
in the one-phonon approximation and from the full 
dynamic structure factor is generally small. Multi- 
phonon processes have clearly the effect of broad- 
ening the spectrum — particularly for large wave 
vectors — but they hardly affect the peaks' posi- 
tions. 

3. In the FM direction at the melting density pi, 
we obtain a substantial improvement over SCP 
resultsji We are thus able to separate the optic 
from the acoustic branches, which appear as dis- 
tinct peaks in the dynamic structure factor. We 
also obtain a significant improvement over previous 
variational QMC rcsultsj^ besides, in Ref. iJ optical 
branches are calculated only in the TA direction. 

4. For the higher density p2 we observe a stronger 
discrepancy between theoretical and experimental 
data, particularly in the FM direction, possibly due 
to the pair potential adopted herci^ 



Transverse modes 

The direct evaluation of transverse phonon modes from 
the dynamic structure factor is hindered by the Q • e(q|j) 
term appearing in its expression (Eqs. [3] and [5]) that 
selects longitudinal modes. At least two strategies can 
be deployed to circumvent this problem. The most im- 
mediate solution consists in considering the peaks in 
the Fourier transform of the transverse counterpart of 
the one-phonon contribution to the dynamic structure 
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Figure 2: (Color online) Longitudinal phonon energies at the melting density p\ extracted from ^(Q,^;) and S'i(Q,(i;), VA and 
VM directions (resp. left and right panel). The limit of the FBZ is indicated by a vertical grey line. Experimental data from 
[27I and [2^. Frequencies calculated at discrete wave- vectors are interpolated by cubic splines as a guide to the eye. 
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Figure 3: (Color online) Longitudinal phonon energies at the density p2 = 0.037 [A ] extracted from S(Q,cj) and S'i(Q,(.tj), 
VA and VM directions (resp. left and right panel). The limit of the FBZ is indicated by a vertical grey line. Experimental 
data from [2§. Frequencies calculated at discrete wave- vectors are interpolated by cubic splines as a guide to the eye. 
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where u\_^ is the transverse component of the atomic 
displacement from equilibrium. Although legitimate in 
principle, this approach is limited to the weak anhar- 
monic regime and only gives access to the positions of 
the peaks, not to their intensities. A better approach, 
which in principle also gives access to peaks intensities, is 
to mimic closely the experimental practice and calculate 
the dynamic structure factor at wave- vectors Q = G 4- q 
such that q is arbitrarily quasi-perpendicular to Q, so 
that a lattice vibration polarized parallel to Q is actually 
quasi-transverses^ this is always possible, just choosing 
a large enough G, i.e. looking at wave- vectors in the 
second, third, or successive Brillouin zones. In such a 
geometry transverse phonon energies at high-symmetry 



wave- vectors in the FBZ can be estimated from the peaks 
in the dynamic structure factor. In order to achieve a 
close comparison between our results and experimental 
data, our simulations at the melting density have been 
performed along the same nuasi-transverse wave-vector 
directions as used in Ref. [271. 

The transverse phonon energies thus obtained are pre- 
sented in Figures [4] and [5l in an extended zone scheme. 

The main findings that emerge from the calculations of 
the transverse modes closely parallel the results obtained 
in the longitudinal case: 

1. The overall agreement with the experimental data 
is good. 

2. The discrepancy between the energy obtained from 
the full dynamic form factor and from its one- 
phonon component is small. Energies from the full 
form factor tend to be more noisy than in the Ion- 
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Figure 4: (Color online) Transverse phonon energies at the melting density pi extracted from the quasi-transverse geometry of 
S{(^,Lu) and from the transverse components of 5'i(Q,c<j), TA direction and Ty branch of the TM direction (resp. left and right 
panel). The limit of the FBZ is indicated by a vertical grey line. Experimental data from[23 and [2^. Frequencies calculated at 
discrete wave- vectors are interpolated by cubic splines as a guide to the eye. 
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Figure 5: (Color online) Transverse phonon energies at the density p2 — 0.037 [A ] extracted from the quasi-transverse 
geometry of S(Q,ti^) and from the transverse components of S'i(Q,a;), TA direction and Ty branch of the TM direction (resp. 
left and right panel) . The limit of the FBZ is indicated by a vertical grey line. Experimental data from [2^. Frequencies 
calculated at discrete wave-vectors are interpolated by cubic splines as a guide to the eye. 



gitudinal case, possibly due to larger multi-phonon 
effects related to the large wave- vector involved in 
the quasi-transverse geometry. 

3. We find a substantial improvement over the SCP re- 
sults, particularly in the TM direction, whilst there 
are no QMC results to compare with. 

4. We find a systematic degradation of the agreement 
with experimental results for increasing density, es- 
pecially for the direction TM of Figure [5l 

B. Intermediate wave-length excitations 

Recent claims that solid *He may display a supersolid 
behavior closely related to superfluidity in the liquid 



phase have prompted a revived interest in the experi- 
mental investigation of density excitations at intermedi- 
ate wave-length, from which valuable information on the 
atomic momentum distribution can be extractedJ^ This 
experimental effort relies on the accurate determination 
of corrections to the Impulse Approximation (Eq. [5]), a 
task which is facilitated in the large-momentum regimefi^ 

Apart from the issues of off-diagonal long-range order 
and Bosc-Einstein condensation in solid "^He, the role of 
atomic interference — as emerging from the additive cor- 
rections to the bare Impulse Approximation — has not yet 
been the subject of detailed theoretical investigations. To 
our knowledge, the best quantitative account of the re- 
sponse of solid helium is limited to the regions of small 
and very large wave- vectors Nonetheless in the inter- 
mediate region where the long- wavelength spectra of both 
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Figure 6: (Color online) Longitudinal excitations along the 
VM direction at the density p\. The solid line is the fit- 
ted free-particle dispersion of a quasi-particle with effective 
mass m*, while the grey line is the free-particle dispersion of 
atomic Helium. The two different background colors ideally 
separate the phonon region from the intermediate region of 
wave- vectors. 



the liquid and the sohd merge into the large-momentum 
regime of nearly free-particle recoil, no ah initio results 
have been reported so far. 
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Figure 7: (Color online) Analysis of the dynamic structure 
factor into the Impulse Approximation and its additive cor- 
rections at the wave-vector Q = 5.92 A along the VM di- 
rection. 

Our QMC methodology, instead, allows us to provide 
an accurate description of this regime as well, showing ev- 
idence of a phonon-like residual coherence, not dissimilar 
to what is found in the superfluid phaseJi We have con- 
centrated our attention to wave-vectors roughly ranging 

from 5 A to 10 A between the phonon and the purely 
single particle regimes. The transition between these two 
regions can be clearly seen looking at the dispersion of 
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Figure 8: (Color online) Coherent oscillations of the cumu- 
lants /i2 and fj.^ appearing in the additive corrections to the 
Impulse Approximation. Cumulants calculated at discrete 
wave-vectors are interpolated by cubic splines as a guide to 
the eye. 

the peaks of the dynamic structure factor, as a function 
of the excitation wave- vector. In Figure [5] the longitudi- 
nal excitation energies of the crystal at the melting den- 
sity pi are shown. In the phonon-like regime (on the left 
of the figure with a yellowish background) we observe 
a periodic dispersion with soft modes corresponding to 
reciprocal-lattice vectors; for larger wave-vectors (on the 
right with a bluish background) we observe a free-particle 
parabolic dispersion, corresponding to an effective mass 
which is slightly higher than the bare Helium mass, being 
M* ~ 1.27 M. 

In order to characterize the density excitations of ^He 
in this region of momenta, we have analyzed the dy- 
namic structure factor in terms of a cumulant expansion, 
Eq. thus extending the scope of previous theoret- 
ical results,— which were essentially unable to predict 
the response of the solid in this intermediate regime. In 
Figure [7] the calculated dynamic structure factor at the 
melting density pi is shown along with its decomposi- 
tion in terms of the leading Impulse Approximation and 
its additive corrections of Eq. PH)) . The coefficients pn 
are taken as free parameters in the fitting procedure of 
the dynamic structure factor, much as it is done in the 
analysis of the experimental datai^^ 

The bare Impulse Approximation of the dynamic struc- 
ture factor overlooks quantum coherence effects in the 
density response functions, which are experimentally 
observediii^ in superfluid '*He. A better account of 
quantum coherence can be achieved including higher or- 
der corrections in the cumulant expansion of Eq. ^ 
whose non- vanishing contribution is in fact recognized in 
the fitting of our reconstructed spectra. Figure [71 The 
presence of quantum coherence between Helium atoms 
in the solid is further appreciated upon looking for devi- 
ations from the incoherent approximation of Eq. A 
cumulant expansion of the incoherent dynamic structure 
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factor can be carried out, and it is known^ that the cumu- 
lants of such an expansion /x™*^ and ^""^ increase mono- 
tonically as whereas /i""^ and /i""^ increase monoton- 
ically as Q^. Characteristic Q— dependent oscillations in 
the ratios and i^l^^^/Q'^ can be therefore exploited 

to infer deviations from the purely incoherent response. 
Such oscillations have been observed in the superfluid 
phasoii and theoretically justified within a T-matrix ap- 
proximation of the He-He atom scatteringfi^ A similar 
behavior has also been recently observed in solid ^Hc,— 
although no satisfactory ab-initio theoretical description 
exists yet. The cumulant dissection of the spectral prop- 
erties extracted from our QMC simulations is a natu- 
ral tool to examine the relics of quantum coherence in 
the intermediate wave- vector region, which closely paral- 
lels the experimental analysis. In Figure [5] wc show the 
Q— dependent oscillations in the ratios and [i^lQ^ 

as found in our analysis of the dynamic structure factor, 
which arc a quite clear manifestation of the residual co- 
herence in the dynamics of solid helium in this intermedi- 
ate region. The quantitative aspects of this analysis may 
be influenced by the quality of the Maximum Entropy 
reconstruction of the spectrum. However, the shift of 
the peak position with respect to the free particle recoil 
frequency should be reliable information, as suggested 
by the good agreement of the calculated and measured 
phonon dispersions previously shown. 

IV. CONCLUDING REMARKS 

In this work wc have demonstrated and successfully 
applied a complete scheme to study the lattice dynam- 
ics of crystalline '''He at zero temperature. Although the 
stochastic nature of QMC methods limits us to explore 
the quantum imaginary time dynamics, we have shown 
that quantitative accuracy can be nonetheless achieved. 



One of the most appealing features of our analysis is the 
possibility to directly parallel the experimental investiga- 
tion based on neutron scattering. The study of the full 
dynamic structure factor has allowed us to describe both 
the phonon and the intermediate wave-length regions of 
excitations. 

At lower density, where the quantum fluctuations sub- 
stantially affect the dynamics, wc have obtained satis- 
factory results for the phonon energies, whereas approx- 
imate quasi-harmonic theories have shown difficulties in 
the accurate determination of the vibrational dispersions. 
An interesting point deserving more research is the effect 
of the adopted pair-potential on the calculated phonon 
branches. Our study suggests a degradation of the agree- 
ment to experimental data at higher density which could 
probably be alleviated upon the inclusion of higher than 
two body terms in the interatomic potential. 

In the intermediate regime of wave-vectors, where the 
density excitations are better understood in terms of cor- 
rections to the small wave-length behaviour, we have 
shown that residual coherence in the quasi-particle exci- 
tations is present as in superfluid helium. Altough much 
of the theoretical and experimental efforts on the merging 
between the phonon and the single particle regimes have 
concentrated on the archetypal quantum solid -^He- we 
believe that an extension of such an analysis to other 
quantum solids such as '^He or molecular hydrogens H2 
would surely be worthwhile. 
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